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Two-component Bose condensates with repulsive interaction are stable when gig2 < gi2 is satisfied. 
By tuning the interactions, we show that the instability corresponding to bose-bose phase separation 
always happens at a higher temperature than corresponding to bose-fermi phase separation happens. 
Moreover, we find both the transition temperature Tdw of supersolid and the coherence peak at ^dw 
are enhanced in the mixtures studied. These will make the observation of supersolid in experiments 
more reachable. 
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1. INTRODUCTION 

Supersolids, a concept simultaneously exhibiting su- 
perfluidity and crystalline order, have been studied in- 
tensely over five decades Theoretically, people 
mainly focus on lattice models of interacting bosons and 
fermions such as the Hubbard model and its various gen- 
eralizations and have obtained many important results by 
numerical analysis [•"), (>]. Experimentally, Kim and Chan 
recently reported they found nonclassical rotational iner- 
tia which should be an direct evidence of supersolid based 
on Lcggett's suggestion in solid "^He [7, s], however, it has 
also been pointed out that this observation may not be 
due to supersolid but due to other reasons, such as an 
increase in shear modulus of bulk solid helium [!)], and 
triggered an intense debate [10, 11]. 

Besides the study of supersolids in condensed mat- 
ter systems, ultracold atoms in optical lattices [12] have 
emerged as a parallel platform witli highly controllabil- 
ity to study supersolids. Trapped Bose-Einstein conden- 
sates with dipole-dipolc interaction can produce a "ro- 
ton" minimum in the excitation spectrum [13-15], and 
this led to the prediction of supersolid upon softening of 
the roton excitation energy [16, 17]. Recently, on the ba- 
sis of ofF-resonant dressing of atomic Bose-Einstein con- 
densates to high- lying Rydberg states, people have found 
the effective atomic interactions resulting from such a Ry- 
dberg dressing can also produce a roton minimum and, 
therefore, provide a clean realization of available model 
for supersolidity [18, 19]. 

In this work, we consider the two kind of bosons are two 
hyperfine state of *^Rb, and the fermions are a hyperfine 
state of ^''K and investigate bose-bose-fermi mixtures in 
a square lattice. For the bose-fermi mixtures subjected 
to a square lattice, it has been pointed out that the den- 



sity wave instability introduced by fermions will establish 
crystalline order, while the condensate bosons exhibit su- 
perfluidity, so a supersolid phase emerges at finite tem- 
perature [17]. For the bose-bose-fermi mixtures studied 
here, besides the density wave instability introduced by 
fermions, there is another instability between the two- 
component bose-condcnsatcs when gig2 = gf^i where gi 2 
are the repulsive intraspecies interaction and 1712 is the 
interspecies interaction [20]. When gig2 > g^2i the bose- 
condensates arc mixed and stable. When 171(72 < <7?2i the 
bose-condensates are unstable and tend to either phase 
separation or collapse depending on gi2 > or < 0. In 
this article, we assume the bose-condensates are initially 
mixed and stable, and we find that bose-bose phase sep- 
aration always happens before bose-fermi phase separa- 
tion when we decrease the temperature. Moreover, we 
find both the transition temperature Tdw of supersolid 
and the coherence peak at fcow are enhanced comparing 
to the bose-fermi mixtures case [17]. 

The article is organized as follows. In Sec. 2, we con- 
sider the two kind of bosons are two hyperfine state of 
^^Rb, and the fermions are a hyperfine state of "^^K and 
investigate bose-bose-fermi mixtures in a square lattice, 
and give the fermionic response in the static limit. In 
Sec. 3, we give the details of the instabilities and different 
phases induced by the instabilities, and give a mean field 
description of the supersolid phase. Some conclusions arc 
obtained in Sec. 4. 
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The Hamiltonian for the bose-bose-fermi mixtures 
takes the form H = Hq + Hint with {a =t, 0, |) 
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dx ^g,^p\'ipl^p^^pf + 52V'|V'|^;i/'i (1) 

where i/'j ^ are the bosonic field operators and i/^g is the 
fermionic field operator. In order to assure the mix- 
tures to be stable, we assume all of the interactions be- 
tween bosons are repulsive, with g^p = Airas^apfi?' /m 
(a, /? = 1, 2. In this work, we use a, (3 to label the inter- 
actions, densities and phases of bosons, and use f, i only 
to label the bosonic operators, moreover, when a = /3 
we only keep a for convenience) and 171(72 > 5^2- 5bf = 
27raBFfi'^/Mj the strength of coupling between bosons and 
fermions, we assume that they are equal for both compo- 
nent of bosons, as, a is the intraspecies scattering length, 
(Is, 12 is the interspecies scattering length, and fi is the rel- 
ative mass. Vq(x) = T4[sin^(7ra;/a) + sin^(7ry/a)] is the 
periodic potential produced by the optical lattice with 
wave-length A = 2a and rria are the mass of bosons and 
fermions. As the bosons are two hyperfine state of ^^Rb, 
it is justified to assume 77i-|- = and Vf ^ Vi for simplic- 
ity in the following. Since the fermions are single compo- 
nent, the interaction between them can be neglected due 
to Pauli exclusion principle. 

In order to obtain the Hamiltonian in momentum 
space, we follow the procedures used in Ref. [17] and ex- 
pand the bosonic and fermionic field operators tpa in the 
forms 



E ^kt,iu;kt,i(a^), 
E CkWk(a;), 



(2) 



keif 



where K denotes the first Brillouin zone, &kt,i and Ck are 
the bosonic and fermionic annihilation operators, while 
Wk-|-,4,(a;) and vi^{x) are the Bloch wave functions corre- 
sponding to a single boson ('f or J,) or fcrmion in the 
periodic potential Va, respectively. Since m-|- ~ mi and 
Vf = Vi, W]g,^{x) should be equal to Wki^x). Therefore, 
we use Wk{x) to denote both of them for convenience. 
Substituting Eq.(2) into Eq.(l) and restricting in the low- 
est Bloch band, we obtain the Hamiltonian in momentum 
space as 



+ E 

{k,k',q,q'} 

N 



{k,k',q,q',a} 



TV 



E "kff^k'o-'-q V 

{k,k',q,q',cr} 



(3) 



where N is the number of unit cells, eF,B^(k) denote 
the energy dispersion of the fermions and bosons, re- 
spectively, while C/bf = .9bf / and C/b.c^ = 
g^f, J dx|{(;|*, with ui(x) and w(x), the Wannier functions 
associated with the Bloch band Wk(x) and Wk(x). In 
a deep optical lattice, the Wannier functions w(x) and 
w (x) are well localized around the minimum of Vq,. As a 
result, the Hamiltonian reduces to a familiar Bose-Fermi- 
Hubbard model, and for ep Bcr(k), only nearest neighbor 
hopping survives. 



eB(q) = 2Jb [2 - cos {q^-a) - cos {qyu)] 
eF(q) = -2Jf [cos (qxa) + cos {qya)] , 



(4) 



where Jb.f is the hopping energy for fermions and bosons, 
respectively. The bosonic dispersion relation implies 
fiB = —4:Jb and the bosons will form a zero-momentum 
Bose-Einstein condensation for sufficiently low tempera- 
ture. The fermionic dispersion relation implies the Fermi 
surface at half-filling Up = 1/2 (where /xp = 0, in this 
work. Tip and ?1b denote the number of particles per unit 
cell) and exhibits perfect nesting for = (7r/a,7r/a) 
and van Hove singularities at k = (0, ±7r/a), (±7r/a, 0). 

Integrating out the fermions produces two effects. To 
the first order in Ubf (in this work, we focus on weak 
interaction, so an expansion in ?7bf and a cut at the 
second order are justified), the fermions simply pro- 
duce a (trivial) shift of the bosonic chemical potential 



Mb, 



Mb a 



[/bf'^f- To the second order in Ubf, the 



fermions provide an effective interaction for the bosons 
which depends on the temperature T of the fermionic 
atom gas. 



{k,k',q,q',(T} 

+ ^ E t^B,.2(T,q-q')6+,fek'T^'qiV4 



{k,k',q,q',} 



with 



C/B„(r,q-q') = [/B, + C/2,x(T,q-q'), 
[/B,i2(r,q-q') = UB.^. + U^BX{T,q-q'). (5) 

The fermionic response in the static limit is given by the 
Lindhard function 



x(T,q) 



dk /[6p(k)]-/[ep(k + q)] 
X Vq ep(k) ~ ep(k + q) +ir]' 



(6) 



where vq = (27r/a)^ is the volume of the first Brillouin 
zone, /(e) = 1/[1 + cxp(e/r)] (^p = at haff filling) 
is just the Dirac-Fermi distribution function. The static 
limit is justified if the fermions are much faster than the 
bosons (Jp >> Jb), as then the fermionic response oc- 
curs on much faster timescales than the movement of 
the bosons, and one can safely neglect retardation effects 



3 



[21]. Using the fermionic dispersion relation Eq.(4), the 
Lindhard function exhibits two logarithmic singularities 
at q = and kow The singularity at q — is purely due 
to the logarithmic van Hove singularity in the density of 
states, and the is due to the combina- 

tion of van Hove singularities and perfect nesting. The 
singularity at q = induces an instability towards a se- 
ries of phase separation, while the singularity at k^w in- 
duces an instability towards density wave formation and 
provides a supersolid phase. The two instabilities are 
competing with each other. 




3. INSTABILITIES AND PHASES 



For the weak interaction, when the temperatures is 
well below the superfluid transition temperature Tkt of 
the bosons, the Lindhard function at q = reduces to 
x{T — > 0, 0) and takes the form [17] 



x(r^o,o) 



deNie)dJ{e) 



No In 



16ci Jp 
T 



(7) 



with N{e) ~ 7Voln|16Jp/e|, iVo = l/(27r2jp) and a = 
2exp(C)/7r « 1.13. As x{T -> 0,0) is always nega- 
tive, the coupling between the bosons and the fermions 
induces an attractive interaction, which is proportional 
to C/pj3x(T', 0), between the bosons (see Eq.(5)). This 
attractive interaction has the effect to reduce the re- 
pulsive interactions UB,ai3 between bosons to C/rff.Q^ = 
f/e.o^ + U^^x{T,0). As a resuh, even [/b,iC/b.2 > U^,^ 
(equivalent to (71(72 > 912) initially, t/rff iC/^f a can be 
tuned to equal to J/^,. by lowering the temperature 
to some value. Moreover, a superfluid condensate at 
low temperatures to be stable requires a positive effec- 
tive interaction C/eff,„ > 0. If we take Ub,i as the en- 
ergy unit, and define the ratios J/b u C^b 12/C^B 1 and 
1 as 7; A and k, respectively. The condition 
^e// 12 ~ ^ef f,iUeff,2 defines the critical temperature 
Tbb.ps for bose-bose phase separation. 



Tb 



16ci Jp exp 



A2 



■7 



7Vok2(1 + 7-2A) 



(8) 



The condition U^tt.a — defines two critical temperatures 



Tb 



and Tb 



for bose-fermi phase separation. 



Tbpi.ps = 16ciJpexp 



16ci Jp exp 



7 



(9) 



When A 7^ 1 and 7, it is directly to show that 
(7 — A^)/(l -f 7 — 2A) is always smaller than min{l,7} 
under the constraint A^ < 7, which is the condition 
that the bose condensates are initially mixed (see Fig.l). 
(7 — A^)/(l-|-7 — 2A) < min{l, 7} indicates when we lower 



FIG. 1: (Color online) A — T Phase diagram. Parameters 
are set as Ib ~ 0.55, A^o^^ — 0.39, 7 — 1.5, Ac ~ y/y- 
For A > Ac, the two-component bosons are phase separation 
initially (BBPS), and supersolid corresponding to bosons-t 
(SSI) emerges when the temperature is below Tijiy. For 
A < Ac, the supersolid phase established by the two Bose- 
condensates appears when the temperature is below Tdw with 
Tbb.ps < Tdw. For A < 1 (1 < A < Ac), the two Bose- 
condensates collapse (separate) when the temperature is be- 
low The, PS with Tbb,ps > Tnw- 



the temperature, the bose-bose mixtures are always eas- 
ier to be unstable and phase separated (or collapse, see 
Fig.l) than the bose-fermi mixtures. Moreover, when A 
gets close to the boundary (7 — A^)/(l -I- 7 — 2A) 
decreases very fast, as a result, Tbb.ps increases exponen- 
tially to values much larger than TOaa;{TBPi,ps, T'bp2,ps} 
and easy to reach in experiments. Therefore, such a bose- 
bose phase separation induced by fermions should be easy 
to be observed in experiments. If we continue to lower 
the temperature after the bose-bose mixtures are phase 
separated, we can expect that bose-fermi phase separa- 
tion will happen and all the components will distribute 
separately in space at last. 

Now, we discuss the second instability induced by 
the singularity in the Lindhard function at k^w. Us- 
ing Eq.(6) and the perfect nesting ep(q-|-kDw)= — eF(q), 
the Lindhard function becomes [17] 



x(r,k 



DWj — 



deN{e) 



tanh(e/2T) 



2 



In 



16ci Jp 



T 



(10) 



The combination of van Hove singularities and perfect 
nesting produces a [InT]^ singular behavior. Such a sin- 
gular behavior can produce a roton minimum at kow. 
Within Bogoliubov theory, the bosonic quasi-particle 
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spectrum becomes 

<±(q) = 4(q)+eB(qK[C/Bi(r,g) + [/B.(r,g)] 
±{eliq)nl [U^,{T,q) ~ UUT,q)f 



4e2(q)n2!73^,(T,<7)} 



(11) 



here we have assumed n^i = = J^b- The induced at- 
traction proportional to [/gpx(r, kow) reduces the energy 
of quasi-particles at kow from a pure-bosonic maximum 
(when L^BF = 0, the maximum of _{q) locates at kow) 
to an induced zero roton minimum (i?B,-(kow)) ~ 0) at 
the critical temperature 



16ci Jf exp 



/t| + 2tB(l+7) + 47-4A2 
2A^oK2(l+7 + tB -2A) 



(12) 



^DWi ^Dw) 



with iB = SJb/^bC/bi- As £'B._(kB 
i?B,-(k = 0) = 0, we can expect the boson modes tikowa 
to become macroscopically occupied just like the boson 
mode below this critical temperature. Comparing 
this result to the one obtained in Ref. [17], 



T 



16ci Jf exp 



V(2 + iB)/AF 



we find Tow is always higher than when param- 
eters appearing in both systems take the same val- 
ues (see Fig.l). Moreover, since Tow depends on 

' ' 2jVol'^(i+7+tB-~2Ar exponentially (12), a small change 
of this term may induce a great change of Tow- There- 
fore, such an enhancement of critical temperature can 
be large. However, Tow can not increase as greatly as 
Tbb.ps, since t^^ has to be larger than a critical value 
^SF-Mi ~ 1/3, below which Mott insulating phase emerges 
and the above picture fails [12]. As a comparison, we cal- 
culate Tow based on the parameters used in Ref. [17] and 
find Tow, max/Tow ~ 1-3 (this ratio goes to the maximum 
when 7 — > 1, in Fig.l, we take 7 = 1.5 just for manifest- 
ing every phase) under the constraint Tow, max > Tbb.ps 
(if bose-bose phase separation happens first. Tow reduces 
to To„, and the enhancement effect of Tow misses). Such 
an enhancement of Tow is of realistic meaning, since the 
lower the temperature is, the harder it is to reach in cold 
atomic experiments. 

For temperatures well below Tkt and Tow, both 
the boson mode fokowa and boa are macroscopically 
occupied, therefore, it is justified to use mean fields 
(^kowa) and (boa) to substitute them. Introducing the 
mean fields (bpg) = V'rWV exp{iipoa) and (fekow«) = 
{Aa /2Ubf) ^/N/rio^ exp{iipQ.) with the constraint nsa = 
+ A'^/{4:noaU^p) and neglecting thermal excitations 
of bosonic quasi-particles [17], we obtain the bosonic den- 
sities as 



Aa COS 9a 
?7bf 



nx Try 

cos — cos 

a a 



(13) 



with 9a = 'Poa ~ fa- The phase difference A9 = 
9i — 02 between the two bosonic density waves determines 
whether they are constructive or destructive. Introduc- 
ing (boa) and (^kowa) to Hamiltonian (3) and neglecting 
terms independent of A^, the Hamiltonian per unit cell 
is given as 



H 

N 



Aj + Al , C/3iA2cos2 0i ^ t/B.A2cos2 02 



2UL 



2UL 



C/B.12A1 A2(cos2 9i + cos2 6*2 + 2 cos(6ii - 6I2)) 



4f/|F 



(14) 



The terms in the first and second lines describe the 
increase in the kinetic and interaction energies of the 
bosons due to the modulation of densities triggered by 
the boson modes fokowai while Hp takes the form 



eF(q) A(0i,02) 
A(0i,02) ep(q') / vv 

(15) 

with a constraint q' = q — kow + K/i (the recipro- 
cal lattice vector K;; ensures the constraint q' G K) 
and A{9i,92) — Aicos^i + A2 cos 6*2. Diagonaliz- 
ing the fermionic Hamiltonian, we obtain the fermionic 
quasi-particle excitation spectrum £'F(k, A) = ±[ep(k) -|- 
(AiCos6'i + A2 cos6'2)^]^/^. To determine the phase 
difference A9, we minimize the thermodynamic poten- 
tial fi(T, Ai, A2, 6*1, (?2) and find a constraint between 
9i and 92'. 9i = 92 = stt, with s an integer. There- 
fore, the phase difference A9 = 0, the two bosonic den- 
sity waves are completely constructive and produce a 
stronger density wave. A stronger density wave makes 
the crystalline order favorable, therefore, such a phase- 
locking effect is favorable to form a supersolid phase. As 
9i = 92 = STT, A+ = Ai + A2 is in fact the gap. In- 
troducing A± = Ai ± A2 and rewriting Eq.(14), the 
self-consistency relations {Oa^^ = 0) take the form 



(1 + 7 + is - 2A)A_ 
1 



(7-l)A+, 



2NoK^ 



(l + ^)+7(l-^) 



2A] 



1 



No Jk ""q 



dk tanh [£;F(k, A+)/2T] 



E^(KA, 



(16) 



Setting A+(Tow) ~ and combining the two equa- 
tions above, we reproduce the critical temperature in 
Eq.(12). This confirms the picture that upon soften- 
ing £'B,_(kow) to zero the bosonic density waves char- 
acterized by (6kowa) 7^ emerge with a breaking of 
the discrete symmetry of the optical lattice be right. 
Furthermore, using the density of states A'^a+ (e) = 



5 



N{^e-^-Al)\e\/^e^ + Al, the gap at T 



becomes 



A+(0) = 32Jpexp 




- 7 + is - 2A) 

(17) 

and the standard BCS relation 2A+(0)/Tow = iir/e^ sa 
3.58 holds. This relation implies that the density wave 
have the characteristic of the superfluid, an evidence of 
supersolid. Therefore, T^w is just the critical tempera- 
ture of supersolid to emerge. 

In experiments, the supersolid can be detected via the 
usual coherence peak of a bosonic condensate in an opti- 
cal lattice. The appearance of a coherence peak at kow 
is a symbol that the supersolid appears. Since the weight 
of this coherence peak is proportional to the number of 
bosons condensed at knw; the larger A_|_ (here equiva- 
lent to (feicDwa)) the sharper the peak is. Therefore, 
based on the similarity of the forms between A-|_(0) and 
^Dw, we find a sharper coherence peak at k^w appears in 
bose-bose-fcrmi mixtures compared to the one appearing 
in bosc-fcrmi mixtures [17] when parameters appearing in 
both systems take the same values. Based on the results 
above, we can make the conclusion that it is more fa- 
vorable to observe the supersolid in bose-bose-fermi mix- 
tures than in bose-fermi mixtures. 



4. CONCLUSIONS 

In this article, we have investigated a bose-bose-fermi 
mixture subjected to a square lattice and found that the 
instability corresponding to bose-bose phase separation 
always happens at a higher temperature than the one 
corresponding to bosc-fermi phase separation. Moreover, 
we find both the transition temperature T^w of super- 
solid and the coherence peak at fcow are enhanced in the 
mixtures studied. These will make the observation of 
supersolid in experiments more reachable. 
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